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1. Introduction 

Thioperoxy (XSO or XOS) radicals exist in a variety of chemical environments, and they have 
as a consequence drawn some interest. HSO, an important species in the chemistry of the troposphere, 
has been examined both experimentally [1-3] and theoretically [4, 5]. The halogenated (X =F, Cl or Br) 
peroxy species and isovalent thioperoxy species have been studied less, but they too are potentially 
interesting because oxidized sulfur species and halogen sources are present in the atmosphere [6]. 
Learning the fate of XSO and XOS radicals is important to understanding the atmospheric oxidation 
chemistry of sulfur compounds. Of these, FSO [7, 8] and C1SO [8, 9] are particularly interesting because 
they have been directly detected spectroscopically. 

Recent studies [1 0, 1 1] in our laboratory on the photochemistry of thionyl halides (X 2 SO; where 
X = F or Cl) have suggested new ways to generate XSO species. The laser-induced photodissociation of 
thionyl fluoride, F 2 SO, at 193 nm and thionyl chloride, Cl 2 SO, at 248 nm is characterized by a radical 
mechanism [1 O, 11], 

X 2 SO -> XSO + X. (1) 

The structure of FSO has been characterized experimentally by Endo el cd. [7] employing microwave 
spectroscopy. Using the unrestricted Flartree-Fock (UHF) self-consistent field (SCF) method, Sakai and 
Morokuma computed the electronic structure of the ground 2 A" and the first excited 2 A' states of FSO 
[12]. Electron correlation was not taken into account in their study. 

In a laser photodissociation experiment, Huber et al. identified C1SO mass spectrometrically [13]. 
C1SO has also been detected in low temperature matrices by EPR [9] and in the gas phase by far IR laser 
magnetic resonance [8]. Although the structure of FSO is known in detail, the only study, experimental 
or theoretical, of C1SO has been an ab initio HFSCF study by Hinchliffe [14], Electron correlation 
corrections were also excluded from this study. 

In order to better understand the isomerization and dissociation dynamics of the radical species, 
we have performed ab initio correlated studies of the potential energy surfaces (PES) of CISO and its 
isomer CIOS at the QCISD^yb-SlG' level of theory [15]. For FSO and FOS, more extensive QCISD/6- 
31 lG(2dl) calculations have been possible, and the results are summarized here. 


2. Structure and spectroscopic properties of FSO and FOS in the ground 2 A" state 

The geometry optimizations and harmonic frequency analyses of the isomers FSO and FOS were 
done with the Gaussian92 suite of programs [16]. The importance of including d- and f-polarization 
functions in basis sets for accurate calculation of the equilibrium geometries of fluorine- and sulfur- 
containing species is documented [17], and adequate sets were employed. The equilibrium geometries and 
harmonic vibrational frequency analyses of the isomer radicals in the ground 2 A" state were determined 
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by the UHF method with QCISD electron correlation correction based on the UHF reference state. A 
“frozen core” exclusion of the inner shells from the QCISD calculations was adopted. The optimum 
geometry of the closed-shell FSCFion was also computed with the same 6-31 lG(2df) basis at the QCISD 
level in order to obtain an accurate estimate of the adiabatic ionization energy of FS O. 

Table 1 shows the optimized geometry of ground state FSO determined at the UHF/6-31lG(2df) 
and QCISD/6-311 G(2df) levels of theory. The results of our calculations on FSO are to be found, 
respectively, in the first and third rows. In their UHFSCF calculations, Sakai and Morokuma[12] 
computed the electronic structures of the ground 2 A" and first excited 2 A' states of FSO. They employed 
4-31 G basis sets and augmented them with s- and p-type Gaussian bond polarization functions. Their 
results are reproduced in the second row of Table 1 for comparison. To our knowledge, there has been 
no ab initio correlated study on this system prior to our own. Endo el al. [7] carried out a microwave 
spectroscopic study of ground state FSO and determined its geometry and fundamental vibrational 
frequencies. The fourth row of Table 1 displays the experimental bond lengths and bond angle. 


Table 1. UHF-QCISD/6-31lG(2df) optimized geometries of ground state FSO and FSO + 



F-S (A) 

S-0 

(A) Z F _ S . 0 (°) 

Total energy (a. u. 

HF This work 

1.564 

1.423 

107.4 

-571.8555 

Sakai 8 

1.560 

1,443 

107.9 


QCISD This work 

1.604 

.454 

108.0 

-572.5015 

Experimcnt b 

1.6023 

.4523 

108.3 


FSO + 





HF This work 

1.472 

1.361 

110.8 

-571.5045 

QCISD This work 

1.510 

1.401 

111.5 

-572.1517 

“Ref. [12] b Experiment: 

Ref. [7] 





The bond lengths and angle of FSO + at the UHF and QCISD levels arc displayed in the last two 
rows of Table 1. Experimental work by Endo et al. [7] as well as ab initio SCF studies have shown FSO 
to be a n-radical with a spin-doublet ground 2 A" state. The unpaired electron is in the SO antibonding 
k* (4a") orbital. FS 0+, absent the unpaired electron in the antibonding orbital, exhibits shortened S-0 
and F-S bonds. 

Accurate estimation of the adiabatic ionization potential of FSO is important to the interpretation 
of mulliphoton ionization experiments. In such experiments, FSO + is produced by ionization via a highly 
excited valence or Rydberg state of FSO, a state which may well have almost the same nuclear 
configuration as does FSO + . At the QCISD level the difference in total energies of FSO (-572.50 15 au.) 
and FSO + (-572.1517 a. u.) gives an accurate estimate (probably within ±0.1 cV with the basis sets 
employed) of the adiabatic ionization potential of FS O. The computed estimate, 9.52 eV, falls into a range 
which indicates that two-photon ionization spectroscopy is feasible with commercially available lasers. 

Table 2 presents the geometry of isomeric FOS radical, determined at the UHF/6-311G(2df) and 
QCISD/6-311G(2df) levels. There are no reported experimental data for the species. With FSO, QCISD 
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correlation correction brings the calculated bond lengths and angle closer to experimental, by 0.04A and 
0.6°, respectively, than they are at the UHF level. In FOS, the effect of electron correlation on geometry 
is greater. A change of O. 15 A in the F-0 bond length is introduced by including correlation. Taking the 
FSO results as a guide, one expects the QCISD geometry of FOS to be accurate to within O.OlA and 1°, 
respectively, in the bond lengths and the angle. The S-0 bond is longer and weaker in FOS than in FSO. 
The local minimum in thetriatomic potential surface which corresponds to FOS lies 83.7 kcal/mole above 
the global, FSO, minimum. In two other thioperoxy radicals (XSO — XOS; where X = H and Cl) which 
have been studied, XSO is also more stable than XOS. In HSO — HOS [5], HSO is the mom stable isomer 
by only 5.4 kcal/mol, whereas ClSO[14]is 42 kcal/mol more stable than CIOS. The uncertainty in the 
reported values of the enthalpy of formation of FSO leads to a range of values for the reaction F( 2 P) + 
SO(X 3 Z ) -» FSO( 2 A"). There is a corresponding range of reaction enthalpies reported, from 75 [18] to 
86 kcal/mol [19], depending on the origin of the value for dissociation of FSO( 2 A")to F( 2 P) and 
SO(X 3 Z~) employed. The product F( 2 P) + SO(X 3 Z~) state is the lowest-dissociation asymptote correlating 
with the 2 A" ground state of FSO or FOS. Assuming the upper limit value to be correct, the local 
minimum corresponding to the FOS isomer lies only a few kcal/mol below the dissociation asymptote 
leading to F( 2 P) + SO(X 3 Z“). 


Table 2. UHF-QCISD/6-3 1 1 G(2df) optimized geometries of ground state FOS. 




F-0 (A) 

S-0 

(A) z F _ s _ 0 O 

Total energy (a. u.) 

HF 

This work 

1.374 

1.589 

110.9 

-571.7135 

QCISD 

This work 

1.526 

1.557 

110.9 

-572.3681 


Table 3 displays the computed harmonic frequencies of FSO and the experimentally observed 
fundamentals [7], The S-0 and F-S harmonic frequencies computed at the QCISD/6-311G(2df) level are 
higher by 20 - 30 cm-’ than the corresponding experimental values. The experimentally observed S-0 
stretching frequency in FSO is larger than the frequency (1 148 cm-') of isolated SO. The S-0 bond length 
in FSO is also shorter by 0.03A than that in diatomic SO (1.481-). The S-0 bond in FSO is 
strengthened by the presence of the electronegative fluorine, which reduces repulsion among the 
nonbonding electrons on oxygen and sulfur. 


Table 3. QCISD/6-311G(2df) harmonic vibrational frequencies (in cm-1) of FSO and FOS. a 


co 


i 


CO, 


CO 


3 


FSO 

Calculated 

1240 

791 

413 


Observed 1 ’ 

1215 

763 

396 

FOS 

Calculated 

911 

491 

362 


3 co l5 co 2 , and co 3 correspond, respectively, to the S-0 and F-S stretching, and the F-S-0 bending 
frequencies of FSO, and, respectively, to the S-0 and F-0 stretches and the F-O-S bend of FOS. 
b Experiment: Ref. 7 
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The FOS harmonic frequencies appear in the last row of Table 3. There is a substantial difference 
in the S-0 stretch of the two isomers; the frequency in FSO is larger by about 300 cm-' than in FOS. 
This difference is consistent with the finding that the S-0 bond in FSO is shorter and stronger than in 
FOS. The computed F-S stretching frequency in FSO is also larger by about 300 cm-’ than the F-0 
stretch in FOS. The low F-0 frequency in FOS indicates that the bond is substantially weaker than the 
F-S bond in FSO and consistent with the fact that the potential energy surface is flat near the local FOS 
minimum region. 


3. Structure and spectroscopic properties of C1SO and CIOS in their ground 2 A" state 

A b initio correlated calculations were performed on the 2 A" Cl-S-0 potential surface in the region 
of the C1SO + CIOS isomerization in order to understand the energetic of the process. About two 
hundred QCISD/6-31 G* and QCISD(T)/6-31 G* calculations were performed with the Gaussian92 system 
to map the surface in the region of interest. 

Analysis of the computed potential surface has yielded minimum energy structures and 
spectroscopic properties of ground state C1SO and CIOS. As with the FSO isomers, the optimum 
geometry of singly ionized, closed-shell C1S0 + was also computed. At the optimum geometry, a single 
QCISD(T) calculation was performed to more accurately estimate the total energy. Tables 4 and 5 show 
the QCISD/6-31G* optimized geometries of ground state C1SO, CISO^ and CIOS. 

C1SO, like FSO, is bent with an unpaired electron in arc* orbital. Similar to the case with FSO, 
removing the unpaired electron in C1SO forms a C1S0 + ion with shorter S-0 and S-Cl bond lengths. 
Mulliken population analysis reveals a rather large formal charge on the sulfur. 


Table 4. QCISD/6-31G* optimized geometries of ground state CISOand C1S0 + . Reported total energies 
are QCISD(T). 



cl-s (A) 

S-0 (A) 

^Cl-S-O ( ) 

Total energy (a. u.) 

C1SO 

2.0864 

1.4957 

109.3 

-932.26341 

ClSO + 

1.9564 

1.4500 

112.0 

-931.92706 


Population analysis reveals some of the differences in charge distribution between CISOand CIOS, 
FSO and FOS, and between the chloro- and fluoro- pairs. In both CISOand FSO there is a relatively 
electropositive atom, S, flanked by a quite negative oxygen and halogen. In FSO the halogen is as 
electronegative as O, whereas in Cl SO it is much less so. Charge flows from the sulfur to the oxygen and 
halogen; in about equal amounts in FSO but more to O than Cl in C1SO. The two molecules are therefore 
somewhat ionic, with large dipole moments. At the HF level the ordering of the dipole moments is FSO 
> C1SO > FOS > CIOS. The S-0 bond order in C1SO is about 7570 greater than the Cl-S. In FSO the 
S-0 bond order is twice that of the F-S bond. 

In CIOS and FOS the bond orders of the two bonds in each are equal, The effect on bonding, 
then, of having the most, rather than the least, electronegative atom in the center of each molecule is to 
reduce the S-0 bond from a double to a single bond. There are differences in the actual details of the 
charge distributions in the two species. In CIOS a very electronegative atom is flanked by two larger 
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atoms of roughly equal electronegativity, while in FOS there is a relatively electropositive atom at one 
end of the molecule, with two highly electronegative atoms bonded to each other at the other end. The 
charge on the sulfur, as the Mulliken gross atomic population, is positive and the same in both molecules. 
It is reduced in magnitude from the values for FSOand C1SO. In FOS, F and O both carry the same 
negative charge, smaller in magnitude than the charges on those atoms in FSO, and reflecting the fact that 
they are bonded to each other. In CIOS the central oxygen carries quite a large negative charge; Cl is 
positive to only a slightly lesser extent than S. 


Table 5. Optimized geometry of CIOS. 



Cl-0 (A) 

o-s (A) 

■^-Cl-O-S 

Total energy (a. u.) 

QCISD 

1.7893 

1.6331 

116.7 

-932.19671* 


* Total energy computed at the QCISD(T)/6-31 G* level of theory. 


Table 6 shows the harmonic frequencies of the two isomers. The S-0 stretch in CIOS is about 
300cm" 1 lower in energy than it is in C1SO, exactly as the corresponding frequencies differ in the 
FSO — FOS pair (see Table 3). The lower frequency reflects the weakness of the S-0 bond in XOS 
relative to XSO which has been cited above. However, lheCl-0 stretch in CIOS is about 100 cm' 1 higher 
in energy than is the Cl-S stretch in C1SO. In this pair of frequencies the C1SO — CIOS pair reverses the 
tendency seen in FSO — FOS. The best explanation of the reversal lies in the contrast in atomic 
electronegativities of the atoms involved in the Cl-S, Cl-O, F-S and F-0 bonds discussed above. 


Table 6. Harmonic frequencies (cm-]) for the C1SO radical and its isomer CIOS determined at the 
QCISD/6-31G* level of theory.* 



CO, 

co 2 

co 3 

C1SO 

1098 

479 

294 

CIOS 

771 

602 

309 


* co,,co 2 , and co 3 correspond, respectively, to the harmonic S-0 and Cl-S stretches and the Cl-S-Obend 
in CISO.For CIOS, they correspond to the S-0 and Cl-0 stretching and Cl-O-S bending frequencies. 
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